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Abstract. A new formulation of statistical mechanics is put forward ac- 
cording to which a random variable characterizing a macroscopic body is 
postulated to be infinitely divisible. It leads to a parametric representation of 
partition function of an arbitrary macroscopic body, a possibility to describe a 
macroscopic body under excitation by a gas of some elementary quasiparticles 
etc. A phase transition is defined as such a state of a macroscopic body 
that its random variable is stable in sense of Levy. From this definition it 
follows by deduction all general properties of phase transitions: existence of the 
renormalization semigroup, the singularity classification for thermodynamic 
functions, the phase transition universality and universality classes. On this 
basis we has also built a 2-parameter scaling theory of phase transitions, a 
thermodynamic function for the Ising model etc. 
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Introduction 

A creation of phase transition theory is one of the most important problem of 
statistical mechanics. While we have reached big success last years we have no still 
solution of the problem [16], [19]. 

In this paper we propose a new approach to this problem which is based on 
analysis of a random variable characterizing a macroscopic body. According to 
traditional exposition of statistical mechanics this random variable is to be a me- 
chanical additive integral of motion. Generalizing this statement we postulate 
that random variables in statistical mechanics must be infinitely divisible. Using a 
canonical representation by Levy-Khintchine for a characteristic function of infinite 
divisible random variable we have got a parametric representation of partition 
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function of any macroscopic body and therefore a possibility to analyze its ther- 
modynamic properties. Since a generalized Poisson random variable, which from 
physical point of view is nothing else than the Boltzmann ideal gas random variable, 
can approximate infinite divisible random variable as good as we wish we can 
explain why it is possible to imagine any excited macroscopic body as a gas of 
some quasiparticles or elementary excitations. A general representation of a set of 
random variables characterizing macroscopic bodies in statistical mechanics follows 
from the fact that a set of infinite divisible random variables is a convex set having 
the normal and Poisson random variables as extreme points. 

Analysis of modern physical representations on a nature of phase transitions 
leads to an idea that the phase transition is such a state of a macroscopic body 
that its random variable is stable in sense of Levy. Taking this as a definition 
we has become able to deduce all well known properties of phase transitions. 
The stability property itself uses in its definition as a matter of fact a notion 
of renormalization semigroup which is well familiar in statistical mechanics. A 
number of finite moments for a stable random variable or, which is the same, a 
number of finite derivatives of a partition function for a macroscopic body, is no 
more than the so called characteristic exponent a, < a < 2, which plays role of 
a critical index for the phase transition. It means that in the nature there exist 
phase transitions only of the first and second order. Using a notion of attraction 
domain for the stable random variable we has been able to explain a meaning of a 
phase transition universality and define universality classes. Distribution functions 
of random variables describing macroscopic bodies in a region of phase transition 
are appeared to be also universal and are expressed in terms of the Fox functions. 

We have applied the above theory to a case when macroscopic body is characte- 
rized by a 2-dimensional random variable (e.g. an energy, an order parameter). We 
have succeed to reproduce all results of a 2-paramctric scaling theory and express 
all critical indexes in terms of appropriate characteristic exponents. In a case when 
one of characteristic exponents is equal to 1 the theory leads to well known results 
of the exactly solvable Ising model with field. 

The paper is organized as follows. In chapter 1 we postulate that the random 
variable characterizing any macroscopic body is infinite divisible and consider dif- 
ferent consequences of this postulate. In chapter 2 we give the definition of phase 
transition and then reproduce all characteristic properties of phase transitions at 
first for an elementary case of the 1-dimensional random variable and further for 
a more general case of the many-dimensional random variable. In chapter 3 we 
consider an example of application of general theory to a particular case of 2- 
dimcnsional random variable and as a result deduce a 2-parametric scaling theory 
of phase transitions. At last we conclude the paper with some remarks about 
general meaning of our results. 

1. Infinitely divisible random variable describes macroscopic body in 

statistical mechanics 

A statistical mechanics studies the properties of macroscopic bodies which con- 
sist, according to definition, of an immense number (e.g. the Avogadro number) of 
some identical particles (e.g. atoms, molecules). As result of that the main laws 
of statistical mechanics are laws of a collective, a crowd and every macroscopic 
body is described by a random variable. The random variable X is certainly to be 
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additive and an appropriate probability measure fj,(x) is surely to be multiplicative 
relative to a convolution *, i.e. if wc divide a macroscopic body into n independent 
subbodies then 

(1.1) X = X 1 + ---+X n , 

(1.2) fi(x) = n{xi) * ■ ■ ■ * fi(x n ). 

The described above scheme of partition and successive construction of a macrosco- 
pic body out of its subparts (we call it the algorithm of Democrites of Abdera) is 
used constantly in statistical mechanics. Of course we suppose that such a partition 
is possible for every positive integer n. It means that a following postulate is valid: 

a random variable corresponding to a macroscopic body is infinitely divisible. 

This postulate generalizes and takes place of a similar statement used in the 
traditional exposition of statistical mechanics due to Gibbs according to which the 
random variable of macroscopic body must be a mechanical additive integral of 
motion. We consider this replacement reasonable since in our opinion we have 
to formulate the general principles of statistical mechanics without any reference 
to other parts of physics (e.g. the classical or quantum mechanics). In essence 
this postulate means that if the atom or molecule dimensions are negligible we 
may apply the Democrites algorithm of decomposition of a macroscopic body into 
subbodies endlessly. 

We can formulate this postulate as 

Definition 1. A macroscopic body in statistical mechanics is the physical structure 
which is described by an infinitely divisible random variable. 

Further we are going to study random variables taking values in M. k . We shall 
consider every random variable to be full i.e. we shall suppose that a support of an 
appropriate probability measure do not belongs to a proper hyperplane of R k . By 
changing an origin of coordinates and a unit of measurement we undergo a random 
variable X to an affine transformation X i— » AX + b, where A is a nondegenerate 
linear transformation and b is a vector. Random variables are said to be of the 
same type if they are connected by an affine transformation. Since this relation is 
reflexive, symmetric and transitive a set of random variables is divided into mutually 
disjoint types. 

A probability measure of a macroscopic body is known to have such a form 
(Gibbs, J.W., 1902) [6] 

(1.3) ii{dx) = Z _1 (u)exp(- < u,x >)D(x)dx, 

where < u,x > is a linear functional on a space of random variables, D(x) is a 
density of states of the macroscopic body for a given value x of a random variable 
and 

(1.4) Z{u) = j cxp(- < u,x >)D{x)dx 

is a partition function. 

A logarithm of a partition function is well known to be a generating function 
for cumulants of the random variable and besides of that defines thermodynamic 
functions of the macroscopic body. If for example the additive random variable is 
an energy then 

(1.5) F(u) = -u~ 1 In Z(u) 
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is the Helmholtz free energy and u" 1 is a temperature. 

Since it is better to describe properties of infinitely divisible random variables in 
the language of characteristic functions we state first of all a connection between 
the characteristic function 



(1.6) fi(y) = M cxp(i(y,x)) = J cxp(i(y,x))n(dx) 

and the partition function Z(u). 
Lemma 1.1. 

(1.7) In fj,(y,u) = ]nZ(u — iy) — lnZ(u). 

Proof. Proof consists of substitution of the expression for macroscopic body pro- 
bability measure (1.3) in the definition of characteristic function (1.6). □ 

It follows from the lemma that the characteristic function p,(y, u) after analytical 
continuation in the variable y to the complex space C fc has such an expression 

(1.8) In p,(iv, a) = lnZ(u + v) -LnZ(u). 

If we do not want to use results of the traditional approach to statistical mecha- 
nics we should consider the lemma and therefore the formula (1.8) as a postulate. 

According to Levy P. and Khintchine A.J. [8] the logarithm of characteristic 
function of an infinitely divisible random variable has such a canonical form 

In fi(y) = i{y, a) - (l/2){y,Ry) + 

(1.9) + / (exp (»<!,,*» - 1 - 7^2) M{dx), 

where a e M fe , R is symmetric non- negative bounded operator in R k and a so called 
Levy measure M satisfies a condition 



/ (1 A \\x\\ 2 )M{dx) < +oo. 



Since the canonical form is unique we may wright down an infinitely divisible 
probability measure in a following way 

(1.10) n=[a,R,M]. 

The canonical form of the logarithm of characteristic function for the infinitely 
divisible random variable is valid under analytical continuation in the variable y 
to any containing the point convex tubular domain G C C k and an appropriate 
integral in this formula converges absolutely and uniformly in any compact K C G. 
Let us remind to this end that the domain G C C fe is named as the tubular one 
with base B C R fc if it is of the form 

G = {z EC k \Sz e B}. 

Now we can prove a following 

Theorem 1.2. A logarithm of the partition function of any macroscopic body has 
such a canonical form 

In Z(u + v) = In Z(u) - (v, a) + (1/2) (v, Rv) + 

(1.11) +/ (ex V (-{v,x))-l + -^^)M(dx), 

J«i=\(n\ \ 1 + \\ x \\ J 
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where a e M. k , R is symmetric non-negative bounded operator in R k and M is a 
Levy measure. 

Proof. It follows from formulae (1.8) and (1.9). □ 

This parametric expression for the logarithm of partition function of macroscopic 
body (and therefore for its thermodynamic functions) must be of extreme impor- 
tance in statistical mechanics. It is obvious that different classes of parameters 
must describe different classes of macroscopic bodies and even different classes of 
microscopic theories that lay in their foundations. Let us remind in this connection 
that the quantum mechanics has appeared for the first time as a new kind of 
distributions for a special case of a black radiation. In next section of the paper 
we use the canonical representation of partition function in order to build up the 
phase transition theory. 

An infinite divisibility of random variable characterizing a macroscopic body 
is tightly connected to a scheme of decomposition of the macroscopic body into 
independent subbodies. Let us consider to this end a sequence of decompositions of 
the macroscopic body into subbodies when at a step n we partition the body into 
k(n) subbodies and denote a random variable of the j-th subbody as X n j, 1 < j < 
k(n). Thus we get a sequence of series of independent random variables X n j, 1 < 
j < k(n). Further we shall suppose these random variables X n j, 1 < j < k(n) to 
be infinitely small, i.e. for any e > such a condition will be valid: 

(1.12) lim sup P{\X nj \>e} = 

n ^°° l<j<k(n) 

Then a limit distribution of a sum 

k(n) 
3 = 1 

exists if and only if the random variable is infinitely divisible. 

Let us discuss now a physical interpretation for some classes of infinitely divisible 
random variables. 

A random variable with normal distribution is well known to be infinitely divisib- 
le. It appears as a limit in the sequences of series of independent random variables 
when and only when for any e > is valid a stronger condition 

k(n) 

(1.13) lim Y P{\X n . 3 \ > e} = 

i=i 

We designate conditions (1-12) and (1-13) as the conditions of weak and strong 
fluctuations appropriately. Using this terminology we may state that the random 
variables with normal distribution should characterize weak fluctuations of the 
thermodynamic quantities in a universal way. This is indeed the case in statistical 
mechanics. 

A random variable with generalized Poisson distribution is also known to be 
infinitely divisible. A probability measure for this case has the following form 

oo 

(1.14) /i = exp (-v(M. k )) v m /n\, 

where v means some single-particle measure. In physics we denote the generalized 
Poisson distribution as the Boltzmann distribution which describes an ideal gas of 
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non-interacting particles of classical mechanics. We should remark here that the 
factor n\ in this formula needs special explanation in a traditional derivation of the 
Boltzmann distribution from the Gibbs one by a necessity to take in account the 
identity of particles or to ensure the additivity of thermodynamic functions. In 
our approach all these properties follows from the infinite divisibility of the Poisson 
distribution. 

Let us consider now some consequences of the infinite divisibility of random 
variables which characterize macroscopic bodies in statistical mechanics. 

First of all we should stress a fundamental importance of the Poisson distributions 
in a theory of infinitely divisible distributions: any infinitely divisible distribution 
is a limit of the Poisson distributions. This remarkable property of the Poisson 
distributions means that a following theorem is valid. 

Theorem 1.3. Any excitation of macroscopic body is possible to present as an ideal 
gas of some elementary excitations or quasiparticles with properties depending on 
a type of the macroscopic body. 

Thus we can give for the first time a natural explanation of the fact well known 
in physics on an existence in macroscopic bodies different kinds of quasiparticles 
such as phonons, plasmons, excitons, polarons, magnons etc. This fact as we have 
seen is just a simple consequence of the infinite divisibility of a random variable 
characterizing a macroscopic body. 

Of course if we need to use the quantum mechanics in order to describe physical 
properties of a macroscopic body we must change the Poisson distribution by the 
Bose and Fermi distributions which are also infinitely divisible as it easy to show. 

According to the canonical representation (1.9) a set of infinitely divisible pro- 
bability measures is a convex hull of the normal and Poisson measures (Johanssen, 
S. 1966) [10]. This mathematical fact is based on a fundamental theorem by Krcin 
M.G., Milman D.P. [14] and has also remarkable physical interpretation: probabi- 
lity measures of macroscopic bodies form a convex hull of the probability measures 
of fluctuations and ideal gases. 

Thus in this section we has given a new formulation of the statistical mechanics 
which is based on the infinite divisibility of random variables that characterize 
macroscopic bodies. In a frame of our approach we do not need to use any notions 
of microscopic theories such as e.g. the hamiltonian, moreover our approach form 
natural external bounds for any possible microscopic theories. We have been able 
to wright down the partition function for any macroscopic body and explain a 
possibility to describe any excitation of it by an ideal gas of some quasiparticles. It 
seems that new formulation is simpler and deeper than usual one and more effective 
especially in studies of general questions of statistical mechanics. We demonstrate 
it in the next section. 

2. Stable random variable describes macroscopic body in a state of 

phase transition. 

We have shown above that a random variable of any macroscopic body is infinite 
divisible. Changing the parameter u £ R fc we can change essentially properties of 
the random variable although it continues to be infinitely divisible all the time. In 
particular this random variable can become stable in sense of Levy, P. (Levy, P. 
1924) [15]. Analyzing the modern physical representations on phase transitions we 
have come to a conclusion that the following definition is valid. 
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Definition 2. The phase transition is such a state of macroscopic body when the 
appropriate random variable is stable. 

From this definition we deduce further all well known basic properties of phase 
transitions. But at first in order to make our exposition as simple as possible and 
concentrate on main ideas of the work we consider the case k = 1, i.e. we shall 
suppose that all random variables under consideration take their values in R. 

. Renormalization semigroup. 

A stability of the random variable characterizing a macroscopic body under phase 
transition means the following: if we split the macroscopic body into n subbodies 
with independent and equally distributed random variables Xk, 1 < k < n then we 
can find always such real numbers a n > and b n that 

n 

(2.15) J2 Xk = anX + bn - 

fe=i 

Here the symbol = means an equality of random variables in their distributions. It 
follows easily from this equality that 

a n = n 1/a , < a < 2. 

A quantity a is usually called the characteristic exponent of a stable random 
variable. If a = then the appropriate random variable X is degenerate i.e. it 
takes with a probability 1 only one value. If a = 2 then the random variable has 
normal distribution. A stable random variable X is called strictly stable if b n = 0. 
If a stable random variable X has the characteristic exponent a ^ 1 then there 
exists such a constant c that the random variable X — c is strictly stable. 

The condition (2.15) which characterizes the stability property of a random va- 
riable is defined only for natural numbers n but we can generalize it easily to all 
real positive numbers t. Being expressed in terms of a probability measure fi the 
stability condition (2.15) takes such a form: 

(2.16) = t 1/a /j,*6(b t ), t>0. 

Here the measure /r means the t-th convolution power of the measure /x, the 
measure t x / a fi means the probability measure for the random variable t x ^ a X, i.e. 
t 1 /"^(i?) = ^(t^^-^E) for any Borel set E, and 6(b) means the Dirac measure in a 
point b. 

According to the definition (2.16) a stable probability measure is invariant under 
a transformation semigroup with elements consisting of a convolution power of the 
measure and a subsequent its affine transformation. We designate this semigroup 
as the renormalization one since it is a proper mathematical generalization of 
that used now widely in a modern theory of phase transitions [16], [19]. Thus 
the stability of a probability measure for a case of phase transition means just 
its invariance with respect to the renormalization semigroup. In other words the 
probability measure under phase transition is a fixed point of the renormalization 
semigroup transformations. We want to emphasize here that we have introduced 
the notion of renormalization semigroup without any reference to the order of phase 
transition and therefore it is valid for any type of transition. And it is very essential 
difference of our renormalization semigroup and that used now in the theory of 
phase transition. 

B. Singularities of thermodynamic functions. Classification of phase transitions. 
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A principal characteristic of a stable probability measure is the characteristic 
exponent a. If a = 2 then the stable measure is normal and therefore it has 
moments of all order. If < a < 2 then the stable measure has finite moments 
only of the order 5 where < S < a. The existence of moments of the order 8 
is known to lead to the existence of derivatives of the order 5 for the logarithm of 
characteristic function what is equivalent according to the lemma proven above to 
the existence of derivatives of the order 5 for the logarithm of partition function or 
the thermodynamic function of the macroscopic body. Thus these arguments proves 
that the thermodynamic functions of macroscopic body under phase transition 
have derivatives only of the order S where < 5 < a and the characteristic 
exponent a plays role of the critical index for the phase transition. Moreover the 
thermodynamic functions may have only the power singularities if the characteristic 
exponent < a < 2, a^l and may have only the logarithmic singularity if a = 1. 
We can formulate these arguments more precise. 

Theorem 2.1. The logarithm of partition function under phase transition has 
singularities of such a type: 
for < a < 2, a ^ 1 



(2.17) In Z(u + v) -In Z(u) = c ia r(-a)v a , v>0 

(2.18) = c 2 aT(-a)\v\ a ; v<0 
for a = 1 

(2.19) lnZ(u + v) - \nZ(u) = Civlnv, v>0 

(2.20) = c 2 |u|ln|u|, v<0. 



Proof. After substitution of the canonical expression for infinitely divisible random 
variable (1.9) in the stability condition (2.16) it is easy to prove that the infinitely 
divisible random variable is stable if and only if its Levy measure M in the canonical 
expression has such a form: 

(2.21) M{x) = -c lX - a , x>0 

(2.22) = c 2 (- X y a , x<0, 

where Ci > 0, c 2 > 0, < a < 2. 

It means that the logarithm of characteristic function of the stable random 
variable is a sum of the Fourier transformation 

F [g}(y) = J g(x)exp(i(x,y))dx 
of the distributions a;^ Q_1 , xl" -1 : 

ln/i(y) - lya + daFlx^- 1 }(y) + c^F^I"- 1 }(y). 
If0<a<2,a^l then 

lnfj,(y) — iya + af(-a) (ci exp(— ina/2)(y + i0) a + c 2 exp(ina/2)(y — i0) a ) . 
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After analytical continuation of this expression from a real line y G M to a 
complex plane z = y + iv 6 C we get 

ln/t(z) = iza + c\aT(— a)(— iz) a , > 0, 

= iza + c 2 aP(— a)(iz) a , > 0. 

If a = 1 then 

, n • • , ( ■ V + i0 \ . , /. !/-t0 

In (j,{y) = iya - cny In -i + c 2 iy\n i 

V exp(l - 7)7 V CX P(! - l), 

where 7 = — r'(l) is the Euler constant. After analytical continuation of this 
expression from a real line to a complex plane we get 



\nu(z) — iza — c\iz\n\ , 3z > 0, 

V exp(l-7)7 

= iza + C2izln\ -), 3z < 0. 

Vexp(l -7)7 

Let us take Sz = ±\v\ in accordance with Sz > 0, 3z < 0. At last using the 
lemma established earlier we prove the statement of the theorem. □ 

According to our theorem the logarithm of partition function is continuous 
function of a parameter v. If < a < 1 then the first derivative of the logarithm 
of partition function has a power singularity o"" 1 . If 1 < a < 2 then the first 
derivative of the logarithm of partition function is continuous but the second one 
has a power singularity v a ~ 2 . In statistical mechanics we call a phase transition 
of the n-th order if the first n — 1 derivatives of a thermodynamic function arc 
finite but the n-th derivative is singular (Ehrenfest, P. 1933) [2]. In accordance with 
this classification of phase transitions a range of characteristic exponents < a < 1 
corresponds to phase transitions of the first order and a range 1 < a < 2 corresponds 
to phase transitions of the second order. If a = 1 then the first derivative of the 
logarithm of partition function has a singularity which is a sum of a jump and a 
logarithmic singularity and the second derivative has a power singularity We 
may designate the last type of phase transitions as the A-point. According to our 
theory there are no any phase transitions of the order more than 2. 

Phase transitions of the first and second order differ substantially not only 
quantitatively by the value of their characteristic exponents a but also qualitatively. 
According to the probability theory a contribution to the sum of random variables 
X\ + . . .+X n appearing in the stability condition (2.15) is due to only one summand 
provided that < a < 1 and is due to all summands equally provided that 
1 < a < 2 (see for example [4]). This mathematical result corresponds to the 
point of view generally accepted in physics according to which a phase transition 
of the first order is the result of local germs appearing by fluctuations and a phase 
transition of the second order takes place at once in a whole volume of macroscopic 
body. 

In general situation we can not predict what order of phase transition should be. 
We have only one exception: if a random variable characterizing the macroscopic 
body is bounded then the phase transition is necessarily of the first order. In fact 
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a stable random variable X with parameters 



a and f3 — 



c\ - c 2 
Cl + c 2 



is bounded from the left (from the right) if and only if (3 = 1([3 = —1) and < a < 1 
(Esseen, C.-G. 1965) [3]. In this connection we should remark that an infinite 
divisible random variable may not be bounded from two sides because otherwise it 
becomes degenerated. 

. Universality of phase transitions. Universality classes. 

The thermodynamic function of nacroscopic body depends in a general state 
according to Theorem 1.2 on real numbers a,R and a Levy spectral measure M 
and in a state of phase transition according to Theorem 2.1 only on four real 
numbers a, C\, C2, a. Therefore when the macroscopic body pass from an arbitrary 
general state to a state of phase transition we get essential contraction of a set 
of parameters necessary to describe this body. As a result of that an immense 
set of macroscopic bodies being characterized by different parameters in a general 
state may be described by identical values of a small number of parameters in a 
state of phase transition. Physicists designate this property of phase transitions as 
universality. 

Now we consider the universality of phase transition in details. Let Xi,X 2 , ■ ■ ■ be 
a sequence of independent random variables with a same distribution function F(x). 
If there exist constants A n ,b ni of a such a type that distribution functions of the 
normalized sums S n = A~ x (Y^k=i Xk — b n ) converge to a distribution function G(x) 
when n — > oo we say that the distribution function G(x) attracts the distribution 
function F(x). Further we shall call a set of all distribution functions F(x) attracted 
to the distribution function G(x) as a domain of attraction of the distribution 
function G(x). According to a probability theory the stable distribution functions 
and only they have domains of attractions. When we interpret a stable random 
variable as a random variable of macroscopic body in a state of phase transition then 
we should interpret an appropriate domain of attraction as a class of universality 
of this phase transition. 

According to a probability theory the normal distribution function attracts the 
largest set of distribution functions while domains of attraction of other stable 
distributions consist of those distribution functions which are in some sense similar 
to these stable distributions e.g. they are expressed by the same formulae as these 
stable distributions where however constants have to be changed by slowly changing 
functions (in sense by Karamata, J. 1930) [12]. More precisely the distribution 
function F(x) belongs to the attraction domain of a stable random variable with a 
characteristic exponent a, < a < 2 if and only if for any constant k > 



l-F(x)+F(-x) 



for 



x — > oo; 



1 - F(kx) + F(-kx) 

F(-x) ci 



for 



x — > oo. 



1 - F(x) c 2 



FOUNDATIONS OF STATISTICAL MECHANICS AND THEORY OF PHASE TRANSITIONS! 1 



Under these conditions we may choose in the following way the constants A n , b n 
characterizing the normalization and centering of a random variable: 

(★) for < a < 1 lim nh(A n ) = 1, b n = 0; 

n— >oc 

(**) for 1 < a < 2 lim nh(A n ) = 1, b n = na, for a = MXk] 

x 



(★★★) for a = 1 lim n/i(A„) = 1, 6„ = nA n / ■ 



dF(a;) 



Here h(x) = 1 — F{x) + F(—x). Comparing expressions for the quantity b n in 
cases < a < 1 and 1 < a < 2 we can give another formulation for a qualitative 
difference between phase transitions of the 1-st and 2-nd order: a mean value of 
the random variable is not changed under phase transition of the 1-st order and is 
changed under phase transition of the 2-nd order. 

For constructing of the phase diagrams it may appear useful a notion of a domain 
of partial attraction. Let X\, X2, ... be a sequence of independent random variables 
with the same distribution function F(x). If there exist such a subsequence of 
natural numbers j and constants Aj , bj , that distribution functions of normalized 
sums Sj = Aj 1 (J2i = i Xk — bj) converge to a distribution function G{x) provided 
that j — > 00 then we say that the distribution function F(x) belongs to a domain 
of partial attraction of the distribution function G{x) . Thus a difference between 
domains of complete and partial attraction rests on a difference between a sequence 
of all natural numbers and its proper subsequence. It appears that every infinitely 
divisible random variable (more precisely, the appropriate type) has a domain of 
partial attraction (Khintchinc, A.J. 1937) [13]. Moreover every random variable 
belongs to a domain of partial attraction of only one type or a countable set of 
types or does not belong to any domain of partial attraction (Doblin, W. 1940, [1]; 
Gnedenko, B.V.1940, [7]). If under changing of external parameters a macroscopic 
body can undergo to different phase transitions then the appropriate random variab- 
le must move from a domain of attraction of one stable random variable to another 
one. This movement may take place only through domains of partial attraction. 

D. Density of a stable statistical distribution under phase transition. Rational 
and irrational critical indexes. 

According to probability theory the stable probability measure has density which 
is analytical function with well known power series expansion. Further for simplicity 
of exposition we shall discuss only strictly stable measure fi which canonical repre- 
sentations is defined by the Levy spectral measure only. Since the Levy measure 
of stable probability measure is power function by the formula (2.21) then except 
of the Fourier transformation of this measure it is reasonable to consider also the 
Mcllin transformation of it 



/>oo 

M[g}{s) = M(s) = / g{x)x s - 1 dx. 
Jo 



It is easy to show that the Mellin transformation of stable probability measure has 
the following expression: 

(2.23) '"■M- t 1|r i:;f"';' 

T(ps - p)T(l + p - ps) 

Here we characterize a stable random measure by parameters a, p instead of the 
parameters a, (3 used earlier. A connection between the parameters (3 and p has 
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such a form: 

(3 = cot(7ra/2) tan(7ra(p - (1/2))) 
Now let us carry out an inverse Mellin transformation 

g(x) = (l/2m) / M(s)x- S ds, 
Jc 

where C means a special contour on a complex plane and compare an expression 
thus obtained with the definition of the Fox -ff-function (Fox, C. 1961) [5] in a form 
of contour integral also 

(oi, A{), ... , (a p ,A p ) 
(fti.-Bi), (b q ,B q ) 



Trrrin 



i r nT^nbj+Bjsm^ii-aj-Ajs) 



(2.24) g(x,a,p)=H. 



li 



~ s ds 

2™ Jc n?= n+ i r K + A 3 S ) ni= m+ i(l - h B jS ) 
As a result we get an expression for stable density in terms of the Fox £f-function 
(Schneider, W. 1986) [22] 

(-l/a,l/a), (-p,p) 
(-1,1), (-p,p) 

Using elementary properties of the Fox iJ-function we can show that the stable 
densities for any values of parameters a and p are solutions of some integro- 
differential equations with the integration and differentiation operators in general 
of the fractional order. 

If the Fox f/-function (2.24) has rational parameters Ai, . . . , A p , Bi, ... , B q then 
it degenerates to the Meyer G-function that satisfies the ordinary differential equa- 
tion with polynomial coefficients of the order max(p, q). This differential equation 
has 2 (if p ^ q) or 3 (if p = q) singular points. Thus if the parameters a,p are 
rational numbers a = M /N, p = L/M, (M, N) = 1, then the stable densities are 
solutions of the usual differential equations of the order max(M — 1, TV — 1) and are 
expressed in terms of the Meyer G-function (Zolotarev, V.M. 1994) [24]: 

g(x,M/N,L/M) - (2 7 r) L -( 1 / 2 )( M+A, )(MiV) 1 / 2 a; - 1 x 



(2.25) G N+L _ 2M+L 



1,1,-1 / x iv 



M M 



l_ N-l 1_ L-l 

AT' • • • ' N 1 L ' • - • ' L 

J_ M-l 1_ L-l 

M ' • - • ' Mi Li' ' ' ' L 



As we see the properties of stable probability measures depend essentially on 
whether critical indexes are rational or irrational. In the modern theory of phase 
transitions the models with rational indexes are well known e.g. the Ising model, 
minimal models in 2-dimensional conformal statistical mechanics etc. 

E. Many- dimensional stable random variables and phase transitions. 

The many-dimensional stable random variables are very important from point of 
view of applications. The Levy spectral measure for the many-dimensional random 
variable has much more complicated form than that for the one-dimensional case 
when as we know it is just a power function. And this is one of the main difficulties 
in the many-dimensional case. 

First of all let us give a short review of stable random variables in K fc [11] [9]. 

A random variable X is called stable if there exist a sequence of independent 
and equally distributed random variables {X n } which after normalization it by 
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the linear operators A n 6 GL(R k ) and centering by vectors b n <E M fc converge by 
distribution to a random variable A, 

n 

(2.26) A n Xj + b n ^X. 

3=1 

A full probability measure (i corresponding to the random variable X is stable 
if and only if there exist an operator B and a vector- valued function b(t) of such a 
type that for all t > 

(2.27) n* =t B n*6(b(t)), 
or which is the same 

(2.28) £*(y) =£(*** y)exp(i<y,&(t))) 
(Sakovich, G.N. 1961,1965[20][21]; Sharpe, M. 1969 [23]) 

In accordance with this result an infinitely divisible measure fj, — [a,R,M] is 
stable if and only if for all t > 

(2.29) [ta,tR,t-M] = [a(t),t B Rt B * ,t B M], 



where a(t) = b(t) + t B a + J 



, . ,,,,,, . t B M{dx). 
l + ||a;|| 2 1 + ||f- B x|| 2 J 

The vector b(t) that appears in this formulae has the following form : 

(2.30) b(t) = t ( v- B dv d, 

Ji/t 

where d is some constant vector in R k . (It is possible that this formula is new, we 
have not succeed to find it in the available literature). 
A set 

(2.31) E(ji) = {Be GL(R k )\ f/ = t B ^ * 8{b{t))} 

we shall call a set of exponents for the stable measure fi under consideration. 
A set 

(2.32) S(fi) = {Ae GL(R k )\ A» = n* 6(a) for some a e M fe } 

we shall call a symmetry group for the measure fi. 

A full stable measure \x has one and only one exponent if its symmetry group is 
finite. 

Let £(i?) be a spectrum of the operator B. We can show that 

(2.33) E(B) C {z G C fe | 5Rz > 1/2}. 

Let /i be a full stable measure with an exponent B. Let us present a minimal 
polynom of the operator B in a form of product of polynomials g\ , 172 of such a type 
that a real part of zeros for the polynom g\ is equal to 1/2 and a real part of zeros 
for the polynom g 2 is strictly greater than 1/2. Then we can present the measure 
fi in a form 

(2.34) (j, = m * H2, 
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where \x\ is the measure of normal random variable and the measure /X2 has not 
any normal component. Moreover 

ker gi(B) = supp Hi, ker g2(B) — lin(supp /i 2 ), and 

R fe = supp hi © lin(supp /i 2 ). 

If a full stable measure /x is not normal then its moments 

a p — J \\x\\ p u(dx) 

satisfy the following inequalities : 

(2.35) a p < 00 forO<p<l/A, 

(2.36) a p = 00 for p > 1/A, where 

(2.37) A = max{Mx\ x G £(£)}. 

If 1 ^ E(-B) then there exist a vector a S R fe which depends on B and is of such 
a type that the measure v = /x * (5(a) is strictly stable i.e. 

Each stable measure in R fc is absolutely continuous with respect to the Lebesgue 
measure on R fc and its density has derivatives of all orders. 

3. Example: 2-parametric scaling theory of phase transitions. 

Now let us analyze a simple example of application for the theory put forward 
above in assumption that the random variable characterizing macroscopic body 
takes its values in R 2 . 

Let the stable random variable describing macroscopic body under phase tran- 
sition has an exponent of the following form 

(3.38) B = ^ ^ , < ai < 2, < a 2 < 2. 

We have to study cases 1 ^ E(-B) and 1 6 E(£?) separately. 
(1) At first let be a x ^ 1, a 2 ^ 1. 

Under this condition we have according to the formula (2.30) 

/fl/on _ 1 fl/a 2 _ 1 

(3-39) b{t) = -d u jd 2 

\ 1/ai — 1 1/012 — 1 

where d\,d 2 are some constants. 

Let us take a measure zx in the form 

(3.40) n = v*8\- - . rfl - ) *J 



l/CKl — 1/ \ 1/012 — 1, 

Then substituting this expression for the measure in (2.27) we can convince our- 
selves that the appropriate random variable is strictly stable i.e. v* — t B v. It 
is worthwhile to point out here that a convolution of a probability measure with 
a degenerate one (i.e. the delta-function) is equivalent to changing a frame to 
reference of the random variable. Thus in accordance with the formula (3.40) 
when we change a frame of reference of the random variable then with assumption 
1 S(_B) we can always make the random variable strictly stable that we shall 
suppose further. 
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According to (2.28) the logarithm of characteristic function fi(yi,y 2 ) of the 
strictly stable random variable satisfy the equation 

(3.41) qiHvi,V2) = »(q 1/ai yi, 9 1/a2 2/ 2 ), q > 0. 

This functional equation has general solution (see for example [20]) 

(3.42) n( yi , y2 ) = (vT+vT) v{y?/y?), 

where v(y) is arbitrary function which in the case under consideration is analytical 
with respect to an appropriate local variable. Let us continue the analytical function 
2/2) m a complex plane and introduce generally accepted in physics notations: 
yi = t is a dimensionless temperature, y 2 = h is a dimensionless field, jj-{yi,y 2 ) = 
$(i, h) is a singular part of thermodynamic function under phase transition. In new 
notations the equation (3.41) is nothing else than the equation of scale invariance 
for singular part of thermodynamic function under phase transition which is well 
known in physics. The thermodynamic function itself according to (3.42) takes 
such a form: 

(3.43) $(t, h) = (t ai + h a2 )V(t ai /h a2 ). 

Now let us define the critical indexes of the heat capacity C, the order parameter 
r] and the susceptibility \ which we can express through the thermodynamic func- 
tion in such a way: 

a 2 $ d<S> <9 2 $ 

(3 - 44) C = -^ ^-oh> x = ~dk?- 

Let us consider separately the cases of weak and strong fields. 
(*) The weak fields, t ai > h" 2 . In this case 

(3.45) $(t, h) = t ai f(h/t ai/a2 , sign t), 

where f(x, ±) are real-analytical functions with the following power expansions: 

(3.46) f(x, +) = 1 + fix 2 + f+x 4 + ... 

(3.47) f{x, -) = 1 + frx + / 2 -z 2 + . . . . 

We define the critical indexes a, (3, 7 in the range of weak fields by such a way: 

(3.48) C-|tr a , 77 - \tf for t < 0, X ~ W~ 7 for t > 0. 
(**) The strong fields, t ai < h a2 . In this case 

(3.49) $(t, h) = h a2 g(t/h a2 / ai ), 

where g(x) is real- analytical function with the following power expansion: 

(3.50) g(x) = 1 + g lX + g 2 x 2 + .... 

We define the critical indexes e, S in the range of strong fields by such a way: 

(3.51) C~h- £ , n-h 1 / 5 . 

Now substituting (3.45) and (3.49) in (3.44) and comparing the results of sub- 
stitution with (3.48) and (3.51) appropriately we get the following expression for 
the critical indexes defined above in terms of the characteristic exponents: 

a = 2 — ai, (i = [a 2 — l)a 1 /a 2 , 7 = (2 — a 2 )ai/a 2 , 

(3.52) e = (2 - ai)a 2 /ai, 6 = l/(a 2 - 1). 
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Analyzing asymptotic expansions for the correlation function of order parameter 

(3.53) G(r) = (v(r)v(0)) ~ r~ a exp(-r/r c (i, ft)), r -> oo, 

where r c (t,h) is the correlation radius, physicists introduce three more critical 
indexes : 

(3.54) r c (i, 0) - r c (0, ft) - h~^, ( = o-d + 2. 

We can easily get indexes v, fi from the equation of scale invariance of the singular 
part of thermodynamic function (3.41) in a form 

(3.55) u d <P(t,h) = ^(tu 1 ^, hu 1 ^), 

where u is a unit of measurement of a linear size of macroscopic body and d is a 
dimensionality of a space containing the body. Substituting in this equation the 
thermodynamic function in the form (3.43) we find that 

(3.56) v = a,\jd, fi = a 2 /d. 

Under assumption r c = oo in the asymptotic expression (3.53) for the correlation 
function of order parameter we can show easily that 

(3.57) ( = d(a 2 -2)/a 2 + d. 

The formulae (3.52), (3.56), (3.57), express the critical indexes of phase transi- 
tions in terms of characteristic exponents and reproduce completely results of so 
called the 2-parametric scaling theory of phase transitions. 

We want to draw attention of the reader to some preferred values of the charac- 
teristic exponents. 

For the characteristic exponent ot\ £ (0, 2] the value ct\ — 2 is singled out. 
Under this value of the characteristic exponent the domain of attraction or, which 
is the same, the universality class of phase transition is the largest possible one and 
therefore we should wait for the value a.\ = 2 in majority of phase transitions. 

According to definition a 2 £ (0, 2]. However if we make in addition a very natural 
assumption that the quantity cr = d — 2 + ( = 2d(a 2 — l)/a 2 which appears in the 
expression for the correlation function (3.53) is non-negative then we should accept 
that a 2 £ (1, 2]. If a 2 — > 1 + then a — > +0 and therefore in a critical range where 
r c = oo and G(r) ~ r _<T the correlations of the order parameter will decrease very 
slowly. This circumstance singles out the value a 2 — 1. (We shall not discuss here 
an interesting case when a 2 £ [0,1]. In this case the correlations of order parameter 
in critical range should increase with distance.) 

Let us apply these considerations to analysis of the critical indexes of phase 
transitions which we have in fact in different cases. 

Let us suppose that d — ► oo as it is indeed in the mean field theory. In this case 
limit values of the critical indexes characterizing the correlation function do not 
depend on the characteristic exponents: 

(3.58) v = 0, n = 0, C = -oo. 
If in addition we assume that 

(3.59) ai =2, a 2 = 4/3, 

then other critical indexes take the so called "classical" values 



(3.60) 



a = 0, (3= 1/2, 7=1, e = 0, 6 = 3. 
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Let us consider d = 3. If in addition we assume that 

(3.61) ai = 2, a 2 = 6/5, 
then the critical indexes take rational values 

(3.62) a = 0, (3= 1/3, 7 = 4/3, e = 0, 5 = 5, v = 2/3, fi = 2/5, C = 0, 

which describe well enough numerous experimental data. However the following 
values of the characteristic exponents lead to better agreement with an experiment 

(3.63) ai = 1, 89 ± 0, 05, a 2 = 1, 210 ± 0, 005, 
(2) Now let be a x = 1, a 2 ^ 1. 

Under this assumption we have according to the formula (2.30), 

t i/a 2 _ i 



(3.64) &(*)= tint*, 

V \ l l a v - 1 

where d\, d 2 are some constants. We may put d 2 = if we change in appropriate 
way a point of reference. Assuming that we can present the probability measure \x 
in the following form 

(3.65) \i = v*\ 
where v is a strictly stable measure and 

(3.66) \{yi,y 2 ) = -idiy x \n\yi\. 

Now similarly to the previous case we get for the logarithm of characteristic 
function the following expression: 

(3.67) A(yi,f/a) = -»diyiln|yi| + (yi + y 2 2 ) P (yi/y 2 2 ), 

where p{z) is an analytic function with respect to some local variable. Continu- 
ing analytically the function fi(yi,y 2 ) in a complex plane and ascribing a special 
physical meaning to its independent variables we obtain from jj,(yi,y 2 ) a singular 
part of the thermodynamic function in a vicinity of phase transition. 

If e.g. we assume that a 2 = 16/15, y\ = t 2 , where t is the dimensionless 
temperature and y 2 — h is dimensionless field then the function fi(yi,y 2 ) is trans- 
formed to a singular part of the thermodynamic function of the exactly solvable 
Ising model 

(3.68) $(t, h) = ct 2 In \t\ + (t 2 + h 16 / 15 )^(t 2 /h w ^ 5 ), 

c is some constant. In order to see that it is enough to consider limit cases 

(3.69) h) ~ t 2 (ci In |i| + c 2 ) + c 3 ht 1/8 , t > h 8/15 

(3.70) $(t, h) ~ c 4 /i 16/15 + t 2 (c 5 In \t\ + c 6 ), t < ^ 8/15 , 

where Cj, i = 1, . . . , 6 are some constants. It is well known classical result (Onsager, 
L. 1944, 1947 [17], [18]). 

In other exactly solvable models of statistical mechanics the characteristic expo- 
nent a 2 is a constant or a rational function of the characteristic exponent ol\. 

We obtain similar results for the macroscopic bodies in state of phase transition 
described by random variables with the following exponents 
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4. Conclusion. 

In the paper we have given a new formulation of statistical mechanics and built 
on this basis a theory of phase transition. According to this approach the random 
variables describing macroscopic bodies in statistical mechanics are infinitely divi- 
sible that is just mathematical expression of an atomic-molecular structure of the 
matter. In a state of phase transition these random variables are stable that reflects 
an invariance of the thermodynamic functions of macroscopic bodies with respect 
to transformations of the renormalization semigroup. 

We should point out some questions which have been studied but have not 
presented here: different examples of computation of critical indexes, discussion of 
limit cases (the mean field theory, the Landau theory, the Lee and Yang theory 
etc.), the non-ergodic systems and spin glasses, the non-equilibrium statistical 
mechanics, an interconnections with microscopic theories. We are going to publish 
them elsewhere. 
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